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Abstract 

In this article, we investigate the structures of the pseudoscalar charmo- 
nium and bottomonium in the framework of the coupled rainbow Schwinger- 
Dyson equation and ladder Bethe-Salpeter equation with the confining ef- 
fective potential (infrared modified flat bottom potential). As the current 
masses are very large, the dressing or renormalization for the c and b quarks 
are tender, however, mass poles in the timelike region are absent. The Eu- 
clidean time fourier transformed quark propagator has no mass poles in the 
timelike region which naturally implements confinement. The Bethe-Salpeter 
wavefunctions for those mesons have the same type (Gaussian type) momen- 
tum dependence and center around zero momentum with spatial extension to 
about q 2 = lGeV 2 which happen to be the energy scale for Chiral symmetry 
breaking, the strong interactions in the infrared region result in bound states. 
The decay constants for those pseudoscalar heavy quarkonia are compatible 
with the values of experimental extractions and theoretical calculations. 
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1 Introduction 



Heavy quarkonium, bound state of the heavy quark and antiquark, characterized 
by at least three widely separated energy scales: the hard scale (the mass m of the 
heavy quarks), the soft scale (the relative momentum of the heavy quark-antiquark 
|p| ) and the ultrasoft scale ( the typical kinetic energy of the heavy quark-antiquark 
E), plays a special role in probing the strong interactions in both the perturbative 
and nonperturbative regions. By definition of the heavy quark, m is large in com- 
parison with the typical hadronic scale Aqcd, the corresponding processes can be 
successfully described in perturbative quantum chromodynamics (QCD) due to the 
asymptotic freedom. However, the lower scales |p| and E, which are responsible 
for the binding, can not be accessible by perturbation theory. The appearance of 
multiscales in the dynamics of the heavy quarkonium makes its quantitative study 
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extremely difficult, the properties of the bound states and their decays can provide 
powerful test for QCD in both the perturbative and nonperturbative regions. 

The physicists propose many original approaches to deal with the long distance 
properties of QCD, such as Chiral perturbation theory [I], heavy quark effective 
theory jS], QCD sum rules lattice QCD perturbative QCD jS], coupled 
Schwinger-Dyson equation (SDE) and Bethe-Salpeter equation (BSE) method P, 
nonrelativistic QCD |7j, potential nonrelativistic QCD [8 , etc. All of those ap- 
proaches have both outstanding advantages and obvious shortcomings in one or 
other ways. The coupled rainbow SDE and ladder BSE have given a lot of suc- 
cessful descriptions of the long distance properties of the low energy QCD and the 
QCD vacuum (for example, Refs. EH HU H2] , for recent reviews one can see 
Refs.|131 HI]). The SDE can naturally embody the dynamical symmetry breaking 
and confinement which are two crucial features of QCD, although they correspond 
to two very different energy scales fSJ E]- On the other hand, the BSE is a con- 
ventional approach in dealing with the two body relativistic bound state problems 
[T7j . From the solutions of the BSE, we can obtain useful information about the 
under-structure of the mesons and obtain powerful tests for the quark theory. How- 
ever, the obviously drawback may be the model dependent kernels for the gluon 
two-point Green's function and the truncations for the coupled divergent SDE and 
BSE series in one or the other wavs|18j. Many analytical and numerical calcula- 
tions indicate that the coupled rainbow SDE and ladder BSE with phenomenolog- 
ical potential models can give model independent results and satisfactory values 
[UJ El Ell HH 1121 ECU El- The usually used effective potential models are confining 
Dirac 5 function potential, Gaussian distribution potential and flat bottom potential 
(FBP) O dlini 113 III] • The FBP is a sum of Yukawa potentials, which not only 
satisfies chiral invariance and fully relativistic covariance, but also suppresses the 
singular point that the Yukawa potential has. It works well in understanding the 
dynamical chiral symmetry breaking, confinement and the QCD vacuum as well as 
the meson structures, such as electromagnetic form factors, radius, decay constants 

[HI 1221 122 • 

During the past two years, the experiments have discovered a number of new 
states, for example, the rf c in exclusive B — > KK s K~tt + decays by Belle [21] , the 
narrow D s j states by Babar, CLEO and Belle (22]) evidence for the + (154O) with 
quantum numbers of K + n [2B], and the X(3872) through decay to 7T + 7T~J/tfj by 
Belle [23 • New experimental results call for interpretations, offer opportunities to 
extend our knowledge about hadron spectrum and challenge our understanding of 
the strong interaction; furthermore, they revitalize the study of heavy quarkonia 
and stimulate a lot of theoretical analysis through the charmonia and bottomonia 
have been thoroughly investigated. 

The decay constants of the pseudoscalar charmonium and bottomonium (j] c and 
776 ) mesons play an important role in modern physics with the assumption of current- 
meson duality. The precise knowledge of the those values f Vc and f Vb will provide 
great improvements in our understanding of various processes convolving the r/ c 
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and rjf, mesons, for example, the precess B — > r/cK, where the mismatches between 
the theoretical and experimental values are large (2H1- The 7] c meson is already 
observed experimentally, the current experimental situation with the r\b meson is 
rather uncertain, yet the discovery of the T]b meson is one of the primary goals of 
the CLEO-c research program [21]; furthermore, the rjj, meson may be observed in 
Run II at the Tevatron through the decay modes into charmed states D*D^^Q. It 
is interesting to combine those successful potential models within the framework of 
coupled SDE and BSE to calculate the decay constants of the pseudoscalar heavy 
quarkonia such as r\ c and r/b . For previous studies about the electroweak decays of 
the pseudoscalar mesons with the SDE and BSE, one can consult Refs. jHJ EH 
[TTJ H2J HS1 EI- In this article, we use an infrared modified flat-bottom potential 
(IMFBP) which takes the advantages of both the Gaussian distribution potential 
and the FBP to calculate the decay constants of those pseudoscalar heavy quarkonia. 

The article is arranged as follows: we introduce the IMFBP in section II; in sec- 
tion III, IV and V, we solve the rainbow SDE and ladder BSE, explore the analyticity 
of the quark propagators, investigate the dynamical dressing and confinement, fi- 
nally obtain the decay constants for those pseudoscalar heavy quarkonia; section VI 
is reserved for conclusion. 



2 Infrared modified Flat Bottom Potential 

The present techniques in QCD calculation can not give satisfactory large r behavior 
for the gluon two-point Green's function to implement the linear potential confine- 
ment mechanism , in practical calculation, the phenomenological effective potential 
models always do the work. As in our previous work [THj, we use a gaussian distri- 
bution function to represent the infrared behavior of the gluon two-point Green's 
function, 

4vrG 1 (A; 2 ) = 3 7 r 2 ^e-^, (1) 

which determines the quark-antiquark interaction through a strength parameter w 
and a range parameter A. This form is inspired by the S function potential (in other 
words the infrared dominated potential) used in Refs.|19^ I2()j. which it approaches 
in the limit A — > 0. For the intermediate momentum, we take the FBP as the best 
approximation and neglect the large momentum contributions from the perturbative 
QCD calculations as the coupling constant at high energy is very small. The FBP 
is a sum of Yukawa potentials which is an analogy to the exchange of a series of 
particles and ghosts with different masses (Euclidean Form), 

n 

where N stands for the minimum value of the masses, p is their mass difference, and 
cij is their relative coupling constant. Due to the particular condition we take for 
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the FBP, there is no divergence in solving the SDE. In its three dimensional form, 
the FBP takes the following form: 



" e -(N+jp)r 

V(r) = -5>— — . (3) 

In order to suppress the singular point at r = 0, we take the following conditions: 

V(0) = constant, 
dV(0) = d 2 V(0) = = d n V(0) = Q 
dr dr 2 dr n 

The a,j can be determined by solve the equations inferred from the flat bottom 
condition Eq.(4). As in previous literature [TS1 I2IJ 1221 E3]- n is set to be 9. The 
phenomenological effective potential (IMFBP) can be taken as 

G(k 2 ) = Gi(k 2 ) + G 2 (k 2 ). (5) 



3 Schwinger-Dyson equation 

The SDE can provide a natural framework for investigating the nonperturbative 
properties of the quark and gluon Green's functions. By studying the evolution 
behavior and analytic structure of the dressed quark propagators, we can obtain 
valuable information about the dynamical dressing phenomenon and confinement. 
In the following, we write down the rainbow SDE for the quark propagator, 

n 



/d k X a X a 
j^r i ^ tJL —S{k)^ u —G^{k-p), (6) 



where 



S \p) = iA(p 2 )^ -p + B(p 2 ) = A(p 2 )[i-f -p + mip 2 )}, (7) 
G» v (k) = (5, u -^)G(k% (8) 

and m C) 6 stands for the current quark mass that explicitly breaks chiral symmetry. 

The full SDE for the quark propagator is a divergent series of coupled nonlinear 
integral equations for the propagators and vertexes, we have to make truncations in 
one or other ways. The rainbow SDE has given a lot of successful descriptions of the 
QCD vacuum and low energy hadron phenomena [HIEIEIEIEIj i n this article, we 
take the rainbow SDE. If we go beyond the rainbow approximation, the bare vertex 
7 M 4^ has to be substituted by the full quark-gluon vertex T^(qqg), which satisfies the 
Slavnov-Tayler identity. In the weak coupling limit, g 2 — > 0, two Feynman diagrams 
contribute to the vertex T a ^(qqg) at one-loop level due to the non-Abelian nature of 
QCD i.e. the self-interaction of gluons jHIj. If we neglect the contributions from 
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the three-gluon vertex T a Aggg) and retain an Abelian version, the vertex T^qqg) 
can be taken as ^T^(qqp), where the vertex T^(qqp) is the quark-photon vertex 
which satisfies the Ward-Takahashi identity. In practical calculation, we can take 
the vertex T^(qqp) to be the Ball-Chiu and Curtis- Pennington vertex [321 EH] so as to 
avoid solving the coupled SDE for the vertex T^(qqp). However, the nonperturbative 
properties of QCD at the low energy region suggest that the SDEs are strongly 
coupled nonlinear integral equations, no theoretical work has ever proven that the 
contributions from the vertex Y a ^(ggg) can be safely neglected due to the complex 
Dirac and tensor structures. The one Feynman diagram contributions version of the 
vertex T a ^{qqg) i.e. neglecting the contributions from the vertex T^ggg) in dressing 
the vertex T a ^(qqg) is inconsistent with the Slavnov-Tayler identity [31 . If we take 
the assumption that the contributions from the vertex T a ^(ggg) are not different 
greatly from the vertex T^qqg), we can multiply the contributions from the vertex 
T^qqg) by some parameters which effectively embody the contributions from the 
vertex Tftggg) [H]. 

In this article, we assume that a Wick rotation to Euclidean variables is allowed, 
and perform a rotation analytically continuing p and k into the Euclidean region. 
The Euclidean rainbow SDE can be projected into two coupled integral equations for 
A(p 2 ) and B(p 2 ). Alternatively, one can derive the SDE from the Euclidean path- 
integral formulation of the theory, thus avoiding possible difficulties in performing 
the Wick rotation . As far as only numerical results are concerned, the two 
procedures are equal. In fact, the analytical structures of quark propagators have 
interesting information about confinement, we will make detailed discussion about 
the c and b quarks propagators respectively in section V. 

4 Bethe-Salpeter equation 

The BSE is a conventional approach in dealing with the two body relativistic bound 
state problems [TZj. The precise knowledge about the quark structures of the mesons 
will result in better understanding of their properties. In the following, we write 
down the ladder BSE for the pseudoscalar quarkonia, 

P P 1 fi-7r P rl A k 

S+ 1 (q + j)x(q,P)SZ 1 (q--) = -^- J ^TmXC*, PhuG^q - k), (9) 

where S(q) is the quark propagator, G til/ (k) is the gluon propagator, P M is the four- 
momentum of the center of mass of the pseudoscalar quarkonia, q^ is the relative 
four-momentum between the quark and antiquark, 7 M is the bare quark-gluon vertex, 
and x(q, P) is the Bethe-Salpeter wavefunction (BSW) of the bound state. 

We can perform the Wick rotation analytically and continue q and k into the 
Euclidean region 2 . In the lowest order approximation, the BSW x{QiP) can be 

2 To avoid possible difficulties in performing the Wick rotation, one can derive the BSE from 
the Euclidean path-integral formulation of the theory. 
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written as 

X (q, P) = 75 [iF?(q, P) + 7 • PP 2 °(g, P)+j-qq- PF'(q, P) + i[ 7 ■ q, j ■ P]P 4 °(g, P)] .(10) 
The ladder BSE can be projected into the following four coupled integral equations, 

POO /*7T 

^F^pj^P) = X) / edk / *n 2 6K(iJ), (H) 

the expressions of the H(i,j) and K{i,j) are cumbersome and neglected here. 

Here we will give some explanations for the expressions of H(i,j) . The P(z, j)'s 
are functions of the quark's Schwinger- Dyson functions (SDF) A(q 2 + P 2 /4 + g • P) 
, P(g 2 + P 2 /4 + g • P), A(q 2 + P 2 /4 - g • P) and P(g 2 + P 2 /4 - g • P). The relative 
four-momentum g is a quantity in the Euclidean spacetime while the center of mass 
four- momentum P must be continued to the Minkowski spacetime i.e. P 2 = — m 2 r , 
this results in the g ■ P varying throughout a complex domain. It is inconvenient to 
solve the SDE at the resulting complex values of the quark momentum, especially 
for the heavy quarks. As the dressing effect is minor, we can expand A and B in 
terms of Taylor series of q ■ P, for example, 

A(q 2 + P 2 /4 + q ■ P) = A(q 2 + P 2 /4) + A(q 2 + P 2 /A)'q ■ P + • • • . 

The other problem is that we can not solve the SDE in the timelike region as the 
two point gluon Green's function can not be exactly inferred from the SU (3) color 
gauge theory even in the low energy spacelike region. In practical calculations, we 
can extrapolate the values of A and B from the spacelike region smoothly to the 
timelike region with suitable polynomial functions. To avoid possible violation with 
confinement in sense of the appearance of pole masses g 2 = — m 2 (g 2 ) in the timelike 
region, we must be care in choosing the polynomial functions [20J. For the r\ c meson, 
the mass is about 3.0Ge^, the extrapolation to the timelike region with the quantity 
— m 2 c /4 can be performed easily, however, the large mass of the rjb meson makes the 
extrapolation into the deep timelike region troublesome. Although the r/b meson has 
not been observed experimentally yet, the theoretical calculations indicate that its 
mass is about 9AGeV j^HI- As the dressed quark propagators comprise the notation 
of constituent quarks by providing a mass m(q 2 ) = B(q 2 )/A(q 2 ), which correspond- 
ing to the dynamical symmetry breaking phenomena for the light quarks. We can 
simplify the calculation greatly and avoid the problems concerning the extrapola- 
tions in solving the BSE by take the following propagator for the c and b quarks, 

S-\q 2 ) =ij-q + M C)b , (12) 

where the M Ct b is the Euclidean constituent quark mass with M 2 b = m 2 b (q 2 ) = q 2 
obtained from the solution of the SDE Eq.(6). 

Finally we write down the normalization condition for the BSW, 

j^ Tr Y^ x ^ p)szl + ^ x ^ p) ^) = 2P ^ (13) 

where % = liX + li- 
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5 Coupled rainbow SDE and ladder BSE and the 
decay constants 

In this section, we explore the coupled equations of the rainbow SDE and ladder 
BSE for the pseudoscalar heavy quarkonia numerically, the final results for the SDFs 
and BSWs can be plotted as functions of the square momentum q 2 . 

In order to demonstrate the confinement of quarks, we have to study the ana- 
lyticity of SDFs for the c and b quarks, and prove that there no mass poles on the 
real timelike q 2 axial. In the following, we take the Fourier transform with respect 
to the Euclidean time T for the scalar part (S s ) of the quark propagator [ol 1X31137] . 



where the 3- vector part of q is set to zero. If S(q) has a mass pole at q 2 = —m 2 (q 2 ) 
in the real timelike region, the Fourier transformed S*(T) would fall off as e~ mT for 
large T or log S* = —mT. 

In our numerical calculations, for small T, the values of S* are positive and 
decrease rapidly to zero and beyond with the increase of T, which are compatible 
with the result (curve tendency with respect to T) from lattice simulations jHH] ; for 
large T, the values of S* are negative, except occasionally a very small fraction of 
positive values. The negative values for S* indicate an explicit violation of the axiom 
of reflection positivity |59] . in other words, the quarks are not physical observable 
i.e. confinement. 

For the c and b quarks, the current masses are very large, the dressing or renor- 
malization is tender and the curves are not steep which in contrast to the dynamical 
chiral symmetry breaking phenomenon for the light quarks, m c (0)/m c ~ 1.5 and 
rrib(0)/rhb — 1.1, however, mass poles in the timelike region are absent. At zero 
momentum, m c (0) = 1937MeV and m&(0) = 5105MeV, while the Euclidean con- 
stituent quark masses M c = 1908Mel/ and M& = 5096MeV, which defined by 
M 2 = m 2 (q 2 ) = q 2 , are compatible with the constituent quark masses in the lit- 
erature. From the plotted BSWs (see Fig.l as an example), we can see that the 
BSWs for pseudoscalar mesons have the same type (Gaussian type) momentum de- 
pendence while the quantitative values are different from each other. Just like the 
lighter qq and qQ pseudoscalar mesons [TH], the gaussian type BSWs center around 
zero momentum with spatial extension to about q 2 = lGeV 2 which happen to be the 
energy scale for Chiral symmetry breaking, the strong interactions in the infrared 
region result in bound states. Finally we obtain the values for the decay constants 
of those pseudoscalar mesons which are defined by 




(14) 



(0|9vy 5 <zk(P)> 




(15) 
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here we use 7r to represent the pseudoscalar mesons 3 , 

f Ve = M9MeV; f Vb = 287MeV, (16) 

which are compatible with the results from the experimental extractions and the- 
oretical calculations, f Vc = 335 ± 75MeV(Exp) [H]; f Vc « AOOMeV(Exp) gl]; 
fr, c = 420 ± 52MeV, f Vb = 705 ± 27 MeV(Theor) g2]; f Vc = 292 ± 25MeV(Theor) 
f Vc w 350MeV(T/ieor) jHJ; /„ c = 300 ± 50MeV(Theor) g3]. In calculation, 
the values of m c and rhb are taken as the current quark masses, m c = 1250Me\^ and 
rh b = 4700MeV; the input parameters for the FBP are N = 1.0A, V(0) = -11. OA, 
p = 5. OA and A = 200MeV, which are determined in study of the qq and qQ pseu- 
doscalar mesons |18j . In this article, the Euclidean constituent quark masses for the 
c and b quarks are taken in solving the BSE as the dressing is tender. We borrow 
some idea from the fact that the simple phenomenological model of Cornell potential 
(Coulomb potential plus linear potential) with constituent quark masses can give 
satisfactory mass spectrum for the heavy quarkonia 4 and take larger values for the 
strength parameter w and range parameter A, i.e. w = 2.2GeV and A = 2.9GeV 2 , 
in the infrared region comparing with the corresponding ones used in Ref.|18j. Fur- 
thermore the masses of the pseudoscalar mesons are taken as input parameters. If 
we take the Euclidean constituent quark masses M c = m c (0) and Mb = rrib(0), 
the decay constants for the r\ c and r\b mesons change slightly, f Vc = 357MeV and 
f Vb = 289MeV. 

6 Conclusion 

In this article, we investigate the under-structures of the pseudoscalar heavy quarko- 
nia rjc and r)b in the framework of the coupled rainbow SDE and ladder BSE with 
the confining effective potential (IMFBP). After we solve the coupled rainbow SDE 
and ladder BSE numerically, we obtain the SDFs and BSWs for the pseudoscalar 
heavy quarkonia i] c and i]b. As the current masses of the c and b quarks are very 
large, the dressing or renormalization for the SDFs is tender and the curves are 
not steep which in contrast to the explicitly dynamical chiral symmetry breaking 
phenomenon for the light quarks, however, mass poles in the timelike region are 
absent. We can simplify the calculation greatly and avoid the problems concerning 
the extrapolations in solving the BSE by making the substitution B(q 2 ) — > M and 
A(q 2 ) — > 1. The BSWs for the pseudoscalar heavy quarkonia have the same type 
(Gaussian type) momentum dependence while the quantitative values are different 
from each other. The gaussian type BSWs center around zero momentum with 
spatial extension to about q 2 = lGeV 2 which happen to be the energy scale for 
Chiral symmetry breaking, the strong interactions in the infrared region result in 
bound states. Our numerical results for the values of the decay constants of the 

3 Here we write down the N c explicitly according to the normalization condition Eq.(13). 
4 For an excellent review of the potential models, one can consult Ref . . 
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pseudoscalar heavy quarkonia are compatible with the corresponding ones obtained 
from the experimental extractions and theoretical calculations. Once the satisfac- 
tory SDFs and BSWs for the pseudoscalar heavy quarkonia are known, we can use 
them to investigate a lot of important quantities involving the B, r\ c and rji, mesons. 
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